Excitation spectrum of vortex lattices in rotating Bose-Einstein condensates 
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Using the coarse grain averaged hydrodynamic approach, we calculate the excitation spectrum of 
vortex lattices sustained in rotating Bose-Einstein condensates. The spectrum gives the frequencies 
of the common-mode longitudinal waves in the hydrodynamic regime, including those of the higher- 
order compressional modes. Reasonable agreement with the measurements taken in a recent JILA 
experiment is found, suggesting that one of the longitudinal modes reported in the experiment is 
likely to be the n = 2, m = mode. 

(N : 



t^J- ' Quantized vortices constitute one of the major topics in the study of Bose-Einstein Condensates (BEC) [1]. Follow- 
ing the demonstration of the nucleation of quantized vortices in BEC [2,3], vortex lattices have since been produced 
1 j , 1 experimentally [4-7], and very recently, both the transverse (Tkachenko) and the logitudinal (hydrodynamic) excita- 
tions of the vortex lattices have been observed at JILA [8]. 

Related systems of vortex lattices in type-II superconductors and in superfluid 4 He have been studied extensively 
in the past, resulting in a large number of important theoretical studies [9-12]. Theoretical investigations of the new 
C$ ■ BEC vortex lattices have also been gaining momentum [13-18]. In particular, Anglin and Crescimanno have provided 
an extension to the Feynman-Tkachenko continuum theory for a dense vortex lattice and thereby studied the long 
wavelength transverse excitations at slow rotation taking full account of the finite geometry [19], while Baym [20] 
has calculated the frequencies of Tkachenko waves at general rotation speeds from Thomas-Fermi to the mean-field 
quantum Hall regimes, with the calculated frequencies in good agreement with the observations of Tkachenko modes 
at JILA. In addition, Cozzini, Chevy and Stringari have studied the hydrodynamic excitations of the surface, scissors 
and Kelvin modes [21]. 

In this paper we calculate the frequencies of the logitudinal excitations of vortex lattices by following the procedure 
used by Svidzinsky and Fetter to study low- lying modes of a single vortex line in BEC [22] . This enables us to write 
the excitation frequencies in terms of the frequencies of a vortex- free condensate, and as a corollary we are able to 
' estimate the frequencies of the higher-order compressional modes as well as the surface excitations. 
\ Although a vortex lattice is significantly more complicated than a single vortex, the concept of diffused vorticity, 
^sO • V x v = 2fi of a vortex lattice provides an effective coarse grain averaged description that greatly simplifies the 
formalism by "smoothing out" the effect of individual vortices [12]. The density of vortex lines per unit surface is 
n„ = 20,/ k where ft is the angular velocity of the containing vessel and k — h/M is the circulation. From the area per 
vortex, n^ 1 = irb 2 , one may define the intervortex separation b — yJh/MQ, which provides, along with the healing 
length, a characteristic length scale for the vortex lattice. The healing length £ = l/^/%^^na, where n is the density 
of the gas and a is the s-wave scattering length, provides an estimate of the size of the vortex core. The functional 
I ■ dependence of n v on fl implies that the vortex lines rotate with the containing vessel; it has indeed been observed in 
superfluid 4 He and in alkali BEC's that vortex lattices undergo self-consistent solid- body rotation. 

A description of the dynamics of a system containing a vortex lattice can therefore be more conveniently given 
by transforming to the frame rotating with angular velocity fl. This renders the Hamiltonian time-independent and 
hence facilitates determination of the metastable (rotating) ground state. The free energy in the rotating frame can 
be written in the following form [17]: 
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where V(r) is the confining potential and g = Airah /M is the interparticle coupling constant with the s-wave 
scattering length a. v s b denotes the velocity of solid body rotation, v s b = fl x r, such that v — v s b defines the velocity 
of the fluid in the rotating frame, and —^Mil 2 r 2 represents an effective repulsive centrifugal potential. 

From the Gross-Pitaevskii Equation (CPE) associated with the free energy, one can write the Thomas-Fermi (TF) 

profile for the vortex lattice as I^tfMI 2 ~ \^TF( r )\ 2u2 ( r ) where l^j^l 2 is the TF density in the absence of a vortex 
lattice 
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provides the periodic vortex lattice structure [17]. In any given unit cell of the vortex lattice, assuming a circular 
Wigner-Seitz cell, u 2 {r) has the local form 
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where £ < r < b [17,23]. In the dense vortex limit £ <C b, i.e. the case of large condensate with many vortices for 
which the concept of diffused vorticity works well, one may approximate (u 2 ) ss 1. We therefore write the TF profile 
as 

l.„ (r) |.- a (1- ^ ( 5 ) 

It is noted that the net effect of a vortex lattice is as if the trapping potential has been modified by that with an 
fi-dependent spring constant. With the self-consistent rotation in place, a 3-dimensional (3D) condensate essentially 
takes on a pancake shape, closely approximating a 2D system [8]; we shall therefore consider 2D BEC in this paper. 
From the usual considerations of normalization, the chemical potential of a vortex lattice in 2D is 



Using the GPE under the zero temperature approximation, i.e. almost all particles are assumed to be in the 
condensate, one can write the corresponding hydrodynamic equations of motion with the substitution = n exp[iS(r)] 
where n is the condensate particle density, n = \^\ 2 , and S is the condensate phase from which the condensate velocity 
is defined: v = TiV 'S(r)/M. The resulting hydrodynamic equations in the TF approximation in which the quantum 
pressure term of the form (fi 2 /2my^n)'V 2 y/n is ignored are the well-known continuity and Euler equations of fluid 
dynamics [24]. When transformed to the rotating frame, the equations take the form [12]: 
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_ +V -(nv) = 0, (7) 
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where, on the right hand side, the terms Mv x 211 and Vi-M£l 2 r 2 provide, respectively, the Coriolis and the centrifugal 
forces experienced by the fluid. These terms vanish in the non-rotating frame. In order to solve the equations, we 
expand the velocity field in the rotating frame, v, and the density n to first order in fluctuations. 

Writing n = uq + n' and v = vo + v' where uq and Vo are the equilibrium values, we linearize these hydrodynamic 
equations to determine the harmonic fluctuations in the particle density n! e~ %UJt and the velocity v' e~ %u)t . Following 
the steps in Ref. [24,25] but for the rotating frame, the equations obeyed by n' and v' are found to be 

iwn' = V • (nov'), (9) 

iuv' = |-Vn' + 2fix v', (10) 
M 

where we have made the continuous vorticity approximation for the velocity field in the rotating frame: 

v «0; V x v' w 0. (11) 

This is equivalent to assuming diffused vorticity for the velocity field in the non-rotating frame, Vxv« 217. Equation 
(11) is consistent with the approximation £ <C b, since the deviation of the velocity field from that of the solid body 
rotation occurs mainly near vortex cores. It should also be noted that, if the number of particles TV is small, the 
contribution of the quantum pressure term may be neglected safely only for the lowest multipoles. 

Equations (9-10) may be written in terms of the velocity components v' r and v'i of v' = v' r f+ v'.<\> as follows: 
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We introduce the velocity potential $' where v' = V$', and assume that, in our coarse grained description, the 
normal modes are proportional to e lm ^ i.e. 

n oc exp(im<j>) ; <& oc exp(im<^). (15) 

Using Eq. (15), and integrating Eq. (13) with respect to r, we get 
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Equating the time derivatives of the velocity potential <&', from Eqs. (17) and (18), one gets 
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which can be re-written as a Cauchy equation for the radial component of $' 
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with a solution $' = ciexp(mlnr) + c 2 cxp(— mlnr), where c\ and C2 are arbitrary constants. The derivative with 
respect to r of <J>' gives the velocity in the radial direction 

9$' m 

—— — — [ci exp(mlnr) — c 2 exp(— mlnr)] , (21) 
or r 

and since the radial velocity of the excitation cannot diverge as r — > 0, c 2 = for m > 0, and c\ — for m < 0. 
Changing notations such that c\ = and c 2 = c^~\ one may write 

— = ±— c (±) exp( m lnr) = ±— (22) 
or r r 

where we have explicitly specified the existence of the positive and negative m values. The equation obeyed by the 
velocity potential $' consistent with the assumptions that we have made is finally: 

iu®' = JLn' T 2in&. (23) 
M ' 

Rearranging Eq. (23) and substituting in Eq. (16), one obtains 
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We now introduce the dimensionless units: u>' = Ut rap , H' = Q/cotrap, r ' = t/Rq. where R 2 ^ = 2/xo/M(cj 2 rap — ft 2 ). 
Making the substitutions in Eq. (24), we get 
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+ lu'(lu'± 2Q!)ri = 0. (25) 



3 



In order to calculate the eigenvalues, the density fluctuation n' needs to be specified. The functional form for the 
excitation of the vortex- free condensate is 7Jq oc r l m lp^' ? re l>°)(l_2r" 2 ) where P,[ a ' b \x) is a Jacobi polynomial of order 



n and the scaled position variable r" — r/Ro, with i?o = y 2^ /Mcj| rap [26,27]. Since n! must tend to n ' as fl' — ► 0, 

i.e. the wave function n' should tend to that for the vortex-free case as the rotational frequency tends to zero, and 
assuming the net effect of a coarse grain averaged vortex lattice on the wave functions to be a modification of the 
spring constant for the trapping potential similar to Eq. (5) for the TF ground state, we shall use as our ansatz that 

n'(r')ocr'\ m \pi lmm (l-2r' 2 ). 

The equation obeyed by ti'q is known to be 
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with the eigenvalues given by (w„„) 2 = \m\ + 2n(n + \m\ + 1) [26,27]. With our ansatz for n'(r'), the expectation 
value of Eq. (25) implies 

J{J ± 20') = [\m\ + 2n(n + \m\ + 1)] (1 - O 2 ), (27) 
which may be solved to give the eigenfrequencies for the vortex lattice of a 2D BEC in the rotating frame: 



<±m = \j [H + 2n(n + |m| + 1)] (1 - fi' 2 ) + O' 2 T O'. (28) 

The solution indicates that energy degeneracy is removed for all m in the rotating frame, with the frequency separation 
of 20'. By calculating the spectrum in the rotating frame, the excitation energy is given relative to the metastable 
ground state of rotating vortex lattice. As to be expected, the spectrum reduces to the well-known result for a vortex- 
free BEC as O' — > 0. One may also proceed in the lab frame for comparison with experiments; the corresponding 
rotational hydrodynamic equations in the lab frame are given by Eqs. (7-8) but with the right hand side of Eq. (8) 
having only the Mv x (V x v) term, and v denoting the velocity field in the lab frame. Following through with the 
calculation one obtains the spectrum in the lab frame: 



= yj [\m\ + 2n(n + \m\ + 1)] (1 - O' 2 ) + O' 2 ± (m - 1)0'. (29) 



This is as expected since the free energy in the rotating frame, F, is related to the energy in the lab frame, H, via 
F = H — L z tt where L z = r x p is the ^-component angular momentum operator. The frequency of the m — 
modes is identical in both the rotating and non- rotating frames, and energy degeneracy is restored for the m = ±1 
modes in the lab frame. It should be noted that, due to the nature of the hydrodynamic calculation, the frequencies 
predicted by Eq. (29) correspond to those of the common-mode logitudinal wave of the vortex lattice rather than the 
transverse Tkachenko waves. It is straightforward to see from Eq. (29) that one can write for the surface (n = 0) 
modes uj' :i _ m — \f\m\ — (jraj — 1)£1' 2 ± (m — 1)0', in agreement with the result of Ref. [21]. In particular, the dipolc 

(or the center of mass mode gives u)' ±1 = 1 for all £1', and the quadrupole mode gives Wg.±2 = V2 — fl /2 ± O' which 
accurately predicts the experimental value of the quadrupole (n = 0,m = ±2) excitations [7]. 

In Ref. [8] , experimental measurement of the frequencies of the common- mode longitudinal waves have been reported 
for the condensate rotating at O' = 0.95. Of the three distinct m = modes that have been observed, the nature of 
the two of the modes have been identified while the third mode could not be classified conclusively. The two modes 
that could be identified were the Tkachenko s-bend mode at 0.072^- and a radial breathing mode at 2.0 ± 0.1 j 1 . 
The third mode, which was presumed likely to be another logitudinal mode, was found to occur at 2.2 ± 0.1^. For 
O' = 0.95, one of the frequencies predicted by Eq. (29) for the n — l,77i — radial breathing mode is uj'-y q — 2.09 
while the frequency for the n = 2, m = mode is predicted to be w' 2 — 2.39. These values differ from the reported 
experimental values by 5 and 7% respectively, or even less when the experimental error bars are taken into account. 
We believe therefore that the third m = mode observed in Ref. [8] is the n = 2, m — compressional mode. 

The dispersion curve, to vs. m, is presented in Fig. 1 for £1' = 0, 0.25, n = 0, 1, 2, and m: m = 0, ±1, ±2, ... ± 10 in 
the lab frame. In addition, it is convenient to define the fractional shift of the squared frequency: 



J 2 - (lo ) 2 _ ±2(m - l)»y [|to| + Mn + jmj + 1)] (1 - Q' 2 ) + W 2 
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The fractional shift as a function of the rotation frequency in the lab frame, 0', for n = 0, 1, 2 and m = 0, ±1, ±2 is 
presented in Fig. 2. 

In summary, we have presented the excitation frequencies of the hydrodynamic shape oscillations in a BEC sup- 
porting a vortex lattice. The methodology we used is similar to that originally applied by Svidzinsky and Fetter [22] 
for the single vortex excitations; our result may therefore be viewed as an application of their result to the case of 
vortex lattices. The most pronounced effect of the presence of a vortex lattice is clearly the fi-dependence of the 
excitation spectrum. 

We note that the hydrodynamic amplitudes of the fluctuation in the particle density and the velocity potential may 
be written in terms of the amplitudes u and v of the Bogoliubov equations, as follows: 

n' = (31) 



*'=2MW^ + ^ <»> 

and the linearized Euler and the continuity equations may be derived from the equivalent GPE [22]. 

The equations derived in this paper arc general; the approximations that we have made in order to model a vortex 
lattice, and also to solve the equations are in the fluid velocity field and the TF density profiles of the ground and 
the excited states. By refining each of these three approximations one may expect to improve on the result obtained. 
For instance, it should be noted that, although the TF density is rotationally symmetric, a real vortex lattice state is 
not strictly an eigenstate of L z . In fact, the only vortex state which is an angular momentum eigenstate is the case 
of a single vortex at the symmetry axis, and if the single vortex is moved off-axis it gives a value of L z which is not 
an integer [28]. It is noted that our analytical result Eq. (29) still holds for non-integer values of m since the Jacobi 

Polynomials Pn a ' b \x) are well defined for all values of a, b > — 1, and a = \m\ and b = in our ansatz for n'. The 
excitation spectra presented in Fig. 1 would therefore be modified accordingly when all such refinements are taken 
into account. As is, our solutions are best applicable to low energy excitations of a large condensate with small £1. 
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